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Prospects of WIMP searches using the annual modulation signature are discussed on statistical grounds, intro- 
ducing sensitivity plots for the WIMP-nucleon scalar cross section. 



1. INTRODUCTION 

The annual modulation effect fil] provides a dis- 
tinctive signature for the identification of a Dark 
Matter signal in the direct searches of WIMPS 
through their elastic scattering off the nuclei of a 
detector. Due to this effect the relative velocity 
between the detector and the WIMP Mawellian 
distribution (assumed at rest in the Galactic rest 
frame) is given by: 



dearth = "sun + «orb sin S COS [u)(t - t )] 



(1) 



where v sun is the Sun's velocity in the galactic 
rest frame, v rb — 30 km sec -1 , sin S ~ 0.51 (6 is 
the angle between the Ecliptic and the Galactic 
plane), lu = 2tt/T, T=l year and to ~ 2 nd june. 

2. EXTRACTING THE MODULATION 
SIGNAL 

Given a set of experimental count rates Nik 
representing the number of events collected in the 
i-th day and k-th energy bin, the mean value of 
N^ (expressed in number of counts per unit of 
detector mass, time and interval of recoil energy) 
is: 

< N ik >= M*fc = (2) 
= [bk + S 0:k + S m , k cosu(ti - t )} ■ Wik (3) 

where the bk represent the average background 
while So,fe and S m ,k are the constant and the 
modulated amplitude of the WIMP signal respec- 
tively. The various parameters of the WIMP 



model are contained in So,k and S m ,k- I n par- 
ticular they depend on the WIMP-nucleus elastic 
cross sections a and the WIMP mass raw- The 
Wik = MATiAEk are the corresponding expo- 
sures, where M is the mass of the detector, AEk 
is the amplitude of the k-th energy-bin, while ATi 
represents the i-th time bin (in the following we 
will assume all AT;= 1 day). For simplicity to 
will be omitted in the following equations. 

The general procedure to compare theory with 
experiment is by making use of the maximum- 
likelihood method. The combined-probability 
function of all the collected Nik, assuming that 
they have a poissonian distribution with mean 
values jj,ik, is given by: 



N lk l 



(4) 



The most probable values of m\\r and a maxi- 
mize L or, equivalently, minimize the function: 



y(mvp,c) = — 21ogX — const 

So.k + S r 



2u-2Y j N lk log [b k 



ik 



(5) 

t COSUjti] 



where y, = ^2 ik l^ik and all the parts not depend- 
ing on mw and a may be absorbed in the constant 
because are irrelevant for the minimization. 
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Once a minimum of the likelihood function has 
been found, a positive result excludes the absence 
of modulation at some confidence level probabil- 
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ity. This can be checked by evaluating the quan- 
tity S 2 — y(a — 0) — y(mw, &)min to test the 
goodness of the null hypothesis. In order to study 
the distribution of S 2 we make use of the asymp- 
totic behaviour: 



6 2 ^ X 2 (a = 0)- X 



2 

rnin 



2 (a,m w ) = ^2 



(S m ,k( m w,cr) - X k ) 2 
Var(X k ) 



J2i N ik cos ujU - N k f3 k 
W k (a k ~(3 2 ) 



(G) 



(7) 



/ WikCOSLOti WikCOS 2 U)ti 

where (3 k = ^ Wk , a k = Wk 

and N k = J2i^ik- In the case of absence of 
a modulation effect numerical simulations show 
that the quantity S 2 belongs asymptotically to a 
X 2 distribution with two degrees of freedom. We 
explain this by the fact that once the cross section 
a is set to zero the likelihood function L no longer 
depends on mw (ah the So and S m functions van- 
ish) and this is equivalent to fixing both the pa- 
rameters of the fit at the same time. In the case of 
presence of a modulation, S 2 has the asymptotic 
distribution of a non central % 2 with one degree 
of freedom and with a mean value given by 



<S 2 >-- 



S m , fc (cr, m w ) 2 AE k 
b k + So, k 



MTa + 2 (8) 



where the same days of data taking have been 
assumed for all the energy bins, and the approx- 
imations J2i Ni k cos 2 uiti ~< N ik > X)cos 2 w£i, 
< Ni k >— Wk(b k + So) have been made. In 
Eq.rtq) we have also defined the factor of merit 
a =lj; J2i cos2 uti (a=l in case of a full period of 
data taking) and the terms depending on the f3 k 
have been neglected. 

Since the degree of overlapping between the dis- 
tributions of 5 2 in the two cases of absence and 
presence of modulation depends on < S 2 >, equa- 
tion (JsJ) allows to estimate the needed exposure 
MTa in order to observe a modulation effect with 
a given probability: for instance, < S 2 >=14.9 
(5.6) corresponds to a 90% (50%) probability to 
see an effect at least at the 95% (90%) C.L. Once 
a required < S 2 > is chosen, a sensitivity plot 
may be obtained by showing the curves of con- 
stant MTa in the plane mw~c. 
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Figure 1. Sensitivity plots in the a^-mw plane for 
a germanium detector. 



4. SENSITIVITY PLOTS AND QUAN- 
TITATIVE DISCUSSION 

In Figures 1-2 we discuss the example of a 
Germanium detector with background b=0.01 
cpd/kg/keV (assumed constant with energy) and 
energy thresholds E t ^=2 keV, values not unreal- 
istic, taken into account the recent performances 
of some Ge detectors. Parameter values used in 
the plots are the local halo mass density p=0.3 
GeV/cm 3 , t>i oc =220 km sec^ 1 ( v\ oc is the mea- 
sured rotational velocity of the Local System at 
the Earth's position), the WIMP r.m.s. velocity 
w ?ms = f v \ oc and "sun ^ («ioc + 12) km sec -1 . 

In Figure 1 the sensitivity plots for < S 2 >=5.6 



is shown in the plane mw — a^ n > , where o-W is 
the WIMP cross section a rescaled to the nucleon 
by adopting a scalar-type interaction. The differ- 
ent curves correspond to values of MTa from 10 
kg- year to 100 kg- year in steps of 10 (from top 
to bottom). The closed contour and the cross in- 
dicate respectively the 2er C.L. region singled out 
by the DAMA modulation search experiment and 
the minimum of the likelihood function found by 
the same authors[||. Note that an exposure of 
10 kg-year of a Ge detector of the above-quoted 
performances would explore almost totally the 
DAMA region. 

In Figure 2 we show, as a function of mw, the 
minimal exposures required for the same germa- 
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Table 1 

Summary of minimal exposures, all in kg • year. Values off (in) parenthesis refer to ni oc =220(170) km sec -1 . 
Eth indicates the energy thresholds expressed in keV, b the background (assumed not dependent on energy) in 
cpd/kg/keV. Exposures are estimated for the WIMP mass range 10< mw <1000 unless specified otherwise. 



Exploration of not DAMA region DAMA region 

excluded regions <5 2 =5.6 <5 2 = 15 



Ge, E th = 2, b=0.1 80(50) 


50(25) 


175(90) 


Ge, E th = 12, b=0.01 25(19)* 


50(25) 


190(95) 


TeO-2, E th = 5, b=0.01 40(25) 


40(20) 


150(80) 


Nal, E th = 2, b=0.1 50 f 


180(100) 


660(355) 



* 45 GeV < m w <110 GeV; ^ w <70(125) GeV. 



nium set-up in order for its sensitivity contour 
to lie below the upper limit on implied by 
the exclusion plot obtained in Ref.||] with pulse 
shape discriminated spectra. Since the present 
uncertainty in wi oc can affect the results in a sig- 
nificant way@|, the different values v\ oc = 170, 220 
and 270 km sec -1 are shown by the dotted, solid 
and dashed curves respectively. In each case the 
values b=0,0.01 and 0.1 are given from bottom to 
top. 

Some examples of minimal exposures for other 
target materials are given in Table 1 for v\ oc = 220 
and 170 km sec -1 . The second column of table 1 
shows the exposures necessary to explore the re- 
gions of the m-w-o^ plane below the exclusion 
plot or Ref. The third and fourth columns 
show the lowest values of MTcv that give a sensi- 
tivity plot encompassing all the 2a DAMA con- 
tour for < 5 2 >=5.6 and 15 respectively. The 
experimental thresholds and resolutions assumed 
in Table 1 are close to those already obtained (or 
foreseeable) in Ge, Nal and Te02 detectors. 

A systematic study of sensitivity plots (not 
shown here for lack of space Q) concludes that 
prospects of modulation searches seem promis- 
ing provided that the WIMP signal is not far 
below present sensitivities, the lowest values of 
explorable <j( n > falling in most cases in the typ- 
ical range of few xlO -10 nbarn. An important 
feature of all the plots is that the sensitivity to 
modulation is generally a decreasing function of 
the WIMP mass, the highest sensitivities corre- 
sponding roughly to the interval 10 GeV < mw < 
130 GeV, and depending in a sensitive way on the 
value of the parameter v\ oc . 
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Figure 2. Minimal exposures MTa for the < 
8 2 >=5.6 calculated for a germanium detector with 
threshold energy E t h=2 keV. 
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